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SPECTROSCOPIC EVIDENCE FOR A SUPERMASSIVE BLACK HOLE IN NGC 4486B 
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ABSTRACT 

The stellar kinematics of the low-luminosity elliptical galaxy NGC 4486B have been measured in seeing 
(7* = 0"22 with the Canada-France-Hawaii Telescope (CFHT) and Subarcsecond Imaging Spectrograph. 
Lauer and collaborators have shown that NGC 4486B is similar to M 31 in having a double nucleus. Here 
we show that it also resembles M31 in its kinematics. Like M31, NGC 4486B rotates fairly rapidly near 



the center (V — 76 ± 7 km s 1 at 0"6) but more slowly farther out (V ~ 20 ± 6 km s 



at r 



4"). 



Also, the velocity dispersion gradient is very steep: a increases from 116 ± 6 km s~ 4 at r = 2"- 6" to 
a = 281 ±11 km s _1 at the center. This is much higher than expected for an elliptical galaxy of absolute 
magnitude Mg ~ —16.8: even more than M31, NGC 4486B is far above the scatter in the Faber- Jackson 
correlation between a and bulge luminosity. Therefore the King core mass-to-light ratio, MjLy — 20, is 
unusually high compared with normal values for old stellar populations {MjLy = 4 ± 1 at Mb — —17). 

We construct simple dynamical models with isotropic velocity dispersions and show that they reproduce 
black hole (BH) masses derived by more detailed methods. We also fit axisymmetric, three-integral 
models. Isotropic models imply that NGC 4486B contains a central dark object, probably a BH, of mass 
M. = 6l 3 x 10 8 M Q . However, anisotropic models fit the data without a BH if the ratio of radial to 
azimuthal dispersions is ~ 2 at r ~ 1". Therefore this is a less strong BH detection than the ones in 
M31, M32, and NGC 3115. A dark mass of 6 x 10 s M© is - 9% of the mass Af bulgc in stars; even if 
M, is somewhat smaller than the isotropic value, M./Mbuige is likely to be unusually large. 

Double nuclei are a puzzle because the dynamical friction timescales for self-gravitating star clusters in 
close orbit around each other are short. Since both M31 and NGC 4486B contain central dark objects, 
our results support models in which the survival of a double nucleus is connected with the presence of a 
BH. For example, they support the Keplerian eccentric disk model due to Tremaine. 

Subject headings: black hole physics - galaxies: individual (NGC 4486B) - galaxies: kinematics and 
dynamics - galaxies: nuclei 



1. INTRODUCTION 

NGC 4486B is a low-luminosity El companion of M 87. 
From the ground, it looks like a normal dwarf elliptical 
similar to but slightly brighter than M32. Faber (1973) 
suggests that it has been tidally truncated by M87; it 
may once have been brighter than its present absolute 
magnitude, Mb — —16.8 (assumed distance = 16 Mpc). 

NGC 4486B was added to the CFHT BH survey (see 
Kormendy 1993; Kormendy & Richstone 1995, hereafter 
KR95, for reviews) because Hubble Space Telescope (HST) 
WFPC1 images showed an elongated, asymmetrical center 
(Lauer et al. 1995). If the double nucleus of M31 (Lauer 
et al. 1993) were moved several times farther away, it 
would look similarly asymmetric. Close double nuclei are 
a surprise because their orbital decay timescales are short. 
At least one promising explanation requires a central BH 
(Tremaine 1995; see §5, below). This makes NGC 4486B 
an interesting target for the BH search. 
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HST WFPC2 images now prove that NGC 4486B has 
a double center (Lauer et al. 1996). Therefore the result 
of the present paper - that it also contains a massive dark 
object (MDO) of mass M, > 10 8 M Q - leaves us with two 
cases of double nuclei found in the presence of MDOs. This 
suggests that the two phenomena are related. 

2. CFHT SIS SPECTROSCOPY 

The observations were obtained with the CFHT and the 
Subarcsecond Imaging Spectrograph (SIS). The slit was 
aligned with the two nuclei. Parameters of the reduced 
spectra are given in Table 1. SIS includes a tip-tilt guider; 
by offsetting the guide probe, we centered the galaxy on the 
slit to 1/4 pixel — O'.'OA. An exposure sequence consisted of 
spectrograph rotation, telescope focus, exposures without 
the grism and slit to center the object, an exposure with 
the slit in place but with no grism, the object spectrum, 
another image through the slit but without the grism, one 
with neither slit nor grism, and a comparison spectrum. 
Between the beginning and the end of a spectral exposure, 
the average drift of the galaxy position was 0"08 ± 0"01 
along the slit and 0"015 ± 0"020 perpendicular to it. 

The seeing was measured on the bracketing direct 
images. Brightness profiles of the galaxy along the slit 
were measured on these and on the spectrum; they agree 
(Fig. 1), so the PSF was the same for the spectrum 
and the direct images. Its Gaussian dispersion radius is 
er« = 0!'28±0"02 (FWHM = 0"66±0!'05) for the average of 
four, 1800 s exposures and cr* = 0"22 for the best exposure. 

After CCD preprocessing and cosmic-ray cleaning, the 
Ca II infrared triplet region of each spectrum was rewritten 
on a log A scale. NGC 4486B is small, so sky spectra from 
the outer parts of the object spectra were subtracted. 
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TABLE 1 
PARAMETERS OF SPECTRA 



Parameter 


Value 


Dates of observations 


1996 April 25 - 27 


Slit position angle 


87° 


Slit length 


168" 


Scale along slit 


0"15 pixel -1 


Slit width 


0!'37 


Wavelength range (1024 pixels) . . 


7580-9405 A 


Reciprocal dispersion 


1.781 A pixel" 1 


Reciprocal dispersion 


63.3 km s -1 pixel -1 


Comparison line FWHM 


2.19 pixel 


Instrumental velocity dispersion 


59 km s _1 


Standard star 


. ri Cyg (K0 III) 



3. KINEMATIC RESULTS 

Velocities and velocity dispersions were calculated with 
Bender's (1990) Fourier correlation quotient program. 
Results are given in Table 2 and Fig. 1. The middle section 
of Table 2 is for the cr* = 0"22 spectrum; the top and 
bottom sections are for the average of all four spectra. 

Since both galaxies have double nuclei, it is interesting 
to compare NGC 4486B with M31. However, note that 
M31 is better resolved by the available spectroscopy than 
is NGC 4486B. In M31, the separation of the nuclei is 
1.8 pc = 0"49 (Lauer et al. 1993), comparable to the 
best spectroscopic resolution (FWHM = 0"64, Kormendy 
& Bender 1996; see KR95) and much larger than the 
spectroscopic scale (0"0864 pixel -1 ). In contrast, the 
nuclei of NGC 4486B are separated by 12 pc = 0"15 (Lauer 
et al. 1996), i. e. one spectroscopic pixel and much less than 
the spectroscopic seeing, FWHM = 0"52 (see Fig. 1). 

NGC 4486B is kinematically similar to M31. Both 
galaxies have rotation curves that increase toward the 
center to apparent maximum velocities that are limited 
by seeing. In M31, V max = 157 ± 4 km s" 1 and in NGC 
4486B, V max = 76±7km s _1 . Outside the spinning center, 
V is only 10-20 km s _1 at r ~ 4" in both galaxies. 
(At still larger radii, bulge rotation is important in M31.) 
Similarly, the velocity dispersion at large radii is normal for 
the luminosity of the bulge (sec Fig. 2) but then increases 
rapidly to remarkably large values at the center. Like V m!ix , 
the apparent maximum velocity dispersion is limited by 
seeing; it is 246 ±8 km s" 1 at <r* = 0"27 in M31 and 
281 ± 11 km s- 1 at cr* = 0"22 in NGC 4486B. 

There are also signs of kinematic asymmetries like those 
in M31. The velocity dispersion is higher on the W side 
of the center ("P2" in Lauer et al. 1996) than on the E 
side. The velocity zero radius is displaced by O'Tl ± 0"03 
from P2 toward PI. And the PI side rotates somewhat 
more rapidly. These effects are marginal at the present 
resolution and should be checked with HST spectroscopy. 

The central a is much larger than wc would expect 
from the Faber- Jackson (1976) correlation between a and 
Mb (Fig. 2). The best stellar-dynamical BH candidates, 
M31 (Dressier & Richstone 1988; Kormendy 1988) and 
NGC 3115 (Kormendy & Richstone 1992; Kormendy et 
al. 1996) stand out: their central velocity dispersions are 
above the scatter for other objects. NGC 4486B is similarly 
above the scatter: a — 281 ± 11 km s" 1 is surprisingly high 
for such a low-luminosity galaxy. Even if NGC 4486B has 
been tidally stripped, a would have been above the scatter 



TABLE 2 

NGC 4486B MAJOR- AXIS KINEMATICS 



r 


V 


e(V) 


a 
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1 1 9 

1 1Z 


94 


-2.96 


18 


14 


140 


15 


-2.21 


40 
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119 
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-1.70 


46 


11 


164 


11 


-1.30 


66 


10 


168 


11 


-1.09 


83 


14 


179 


17 


-0.94 


73 


9 


172 


11 


-0.79 
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24 


0.00 
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0.15 


-47 


18 


284 


19 


0.30 
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13 


180 


15 


0.53 


-93 


16 


194 


18 


0.75 


-100 


17 


218 


15 


1.11 


-63 


14 


156 


17 


0.86 


-86 


12 


172 


14 


1.01 


-70 


15 


174 


16 


1.23 


-64 


11 


168 


11 


1.60 


-57 
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144 


11 


2.57 


-40 


10 


113 


12 


3.70 


-10 


11 


110 


10 


5.51 


-26 


21 


98 


35 



of other objects for any plausible progenitor luminosity. Of 
course, the BH galaxies in Fig. 2 have been measured with 
better resolution than most other galaxies. Our comments 
are not meant to imply that these others do not contain 
BHs; they may move upward too when they are observed 
at high resolution. But the scatter of points in Fig. 2 shows 
the range of a values for which the central mass-to-light 
ratio is normal for an old stellar population. With respect 
to these, BH candidates stand out as having large a. 

For a resolved core, the King core mass-to-light ratio 
(Richstone & Tremaine 1986) provides a good indication of 
whether there is extra dark matter near the center. Where 
So is the central surface brightness, r c is the core radius, 
and G is the gravitational constant, 



Kormendy (1993) emphasizes that Sor c ~ constant in 
different galaxies, so it is usually not r c or S that 
determines whether an object has a high M/L ratio. 
Instead, M/L is high or low if a 2 is high or low compared 
to other galaxies. For a = 281 km s _1 , r c ~ 0"23, and 
Tiq corresponding to 14.26 V mag arcsec -2 , the core mass- 
to-light ratio is M/Ly ~ 20. This is much higher than 
values (4 ± 1 at Mb — —16.8) normally observed in old 
stellar populations. Therefore NGC 4486B's high velocity 
dispersion suggests that it contains an MDO. 
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FIG. 1. — Surface brightnesses and kinematics along the 
major axis of NGC 4486B. In the top panel, the lower solid 
line is for the deconvolved HST image (Lauer el al. 1996) 
convolved with our PSF. Filled circles show the brightness 
profile measured from our best spectrum. 



2.6 

in 2.4 
E 

CS 2.2 
b 

rji 2.0 h 



•N4486B 



-| i i i |- 
•N3115 



•■c3-.., : - >•"•/:•■ . 



-20 



-22 



FIG. 2.— Correlation of log cr with M B for 594, E and SO 
galaxies from McElroy (1995) and for faint ellipticals from 
Bender, Burstein, & Faber (1992). Distances are based on a 
Hubble constant of Hq = 75 km s _1 Mpc -1 as in McElroy 
(1995). Filled circles show apparent central velocity dispersions 
for M31 (CFHT: cr* = 0"27), NGC 3115 (HST, aperture = 
0"21), and NGC 4486B (present data). Open symbols are bulge 
dispersions outside the region affected by the MDO. 



4. SPHERICAL ISOTROPIC DYNAMICAL MODELS 

Detailed mass modeling of NGC 4486B is postponed 
until HST spectroscopy becomes available. Here we take 
a first look at whether NGC 4486B contains an MDO by 
constructing spherical, isotropic dynamical models. The 
gravitational potential due to the stars is calculated from 
the brightness distribution, a central MDO is added, and 
the parameters of both are varied until a match is obtained 



to (V 2 



r 2U/2 



Output are the allowed range of stellar 



mass-to-light ratios M/Ly and MDO masses M,. 

The models are constructed as follows. First, we find the 
galaxy's intrinsic light distribution v(r) by deprojecting 
its surface brightness profile. We assume that M/Ly is 
independent of radius except for M., and we derive the 
model potential $(r) by solving the Poisson equation. For 
each assumed M,, we then integrate the Jeans equation, 

to find the intrinsic second-order moment va 2 (r). We 
project this and the zeroth-order moment v(r) along 
each line of sight and convolve the results with the 
spectroscopic seeing. Finally, we average the seeing- 
convolved zeroth- and second-order moments over the bins 
used in the spectroscopy. Dividing each second-order 
moment by the corresponding zeroth-order moment and 
taking the square root yields the predicted dispersion 
profile to be compared with the data. Error estimates 
on M m include only the uncertainty in the a data; they 
do not take into account differences between velocity 
measurements and velocity moments, errors in the 
assumptions (e.g., anisotropy effects), or distance errors. 
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FIG. 3. — Comparison of BH masses obtained using the 
present, simple modeling technique with masses (KR95) derived 
via more detailed models (e.g., including rotation). 

Figure 3 compares BH masses derived as above with 
those obtained from more detailed models (see KR95 for 
references). The highest- resolution ground-based results 
are used; these match the resolution of the present 
data. Only M32 has been updated since KR95 (Bender 
et al. 1996). Figure 3 shows that for most objects the 
present technique reproduces the published BH masses to 
better than a factor of two. 
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5. A SUPERMASSIVE BH IN NGC 4486B 

The machinery of § 4 was used to fit the HST WFPC2 
photometry in Lauer et al. (1996) and the kinematic data 
in Fig. 1. The results are shown in Fig. 4. These 
models imply that NGC 4486B contains an MDO with 
M. ~ 8.7 x 10 8 M©. The corresponding stellar mass-to- 
light ratio is M/Ly — 5.0. As in other MDO detections, 
the most probable interpretation is a supcrmassive BH. 
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FIG. 4. — Fits of spherical, isotropic models to the kinematics 
of NGC 4486B. The solid line is the best-fitting model (M. = 
8.7 X 10 8 M ; M/Ly = 5.0). The dashed lines are "error 
bar" models with M. = 10.5 X 10 8 M ; M/Ly = 3.8 and 
M. = 6.9 X 10 8 M ; M/Ly = 6.3. The dotted line fitted at 
large radii is a no-BH model with M/Ly = 6.4. 

However, this is a weaker BH detection than the best 
stellar-dynamical cases (the Galaxy, M31, M32, and NGC 
3115). The reason is that axisymmetric, three-integral 
maximum entropy models (Richstone et al. 1996) can fit 
the data without a BH provided that the ratio of radial 
to azimuthal velocity dispersions is a r /ag ~ <J r /<J<p — 2 
at r ~ 1". Since rotation is relatively unimportant in 
NGC 4486B, some anisotropy would not be surprising. 
On the other hand, anisotropic models may be unstable 
(Merritt 1987). Also, isotropic model measurements of M. 
have turned out to be close to correct when anisotropic 
models were constructed (see KR95 for a review). E.g., 
in M87, the Sargent et al. (1978) measurement of M, has 
turned out to be remarkably accurate (Harms et al. 1994) 
even though the galaxy is so luminous that anisotropy is 
almost assured and even though anisotropic models can 
explain the kinematics without a BH (Duncan & Wheeler 
1980; Binney & Mamon 1982; Richstone & Tremaine 1985; 
Dressier & Richstone 1990; van der Marel 1994). Therefore 
it is likely that M, in NGC 4486B is close to but somewhat 
smaller than the value given above. Even nearly isotropic 
maximum entropy models give a smaller mass than those of 
Fig. 4, i. e., M. ~ 4x 10 8 M©. In this case, the stellar mass- 
to- light ratio Mj Ly ~ 7 is unusually high for a galaxy of 
Mb = —16.8, while the models of Fig. 4 imply a more 
nearly normal mass-to-light ratio. Therefore: 



We adopt the mean of the Jeans and maximum entropy 
model results, M. = 6±^ x 10 8 M© and M/Ly = 6 ± 1. 

This is a remarkably large BH mass for such a low- 
luminosity galaxy. It is well above the correlation between 
M, and bulge luminosity that is observed for other galaxies 
(Fig. 5). For ten previous BH detections, the mean 
ratio of M, to the mass of the bulge is <M./Mbui gc >= 
0.0022±g;ggJ|. To date, NGC 3115 has had the largest 
BH mass fraction, 2.4%. But 6 x 10 8 M© is 9A± 5 3 ° 5 % of 
the mass in stars. Even if the galaxy once was brighter, 
its BH mass is surprisingly high. In fact, NGC 4486B 
is a good model of what a nearly naked quasar (Bahcall 
et al. 1994, 1995, 1996) would look like after it turned 
off: a monstrous BH in a tiny galaxy. However, since 
anisotropic models allow smaller BH masses, it will be 
important to check this result when HST spectroscopy 
becomes available. 
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FIG. 5. — BH mass as a function of bulge absolute magnitude, 
from KR95 but with BH masses updated for M 32, NGC 3115, 
and NGC 4594, and with NGC 205 (Jones et al. 1996), NGC 
4261 (Ferrarese et al. 1996), NGC 3379 (Gebhardt et al. 1996) 
and NGC 4486B (this paper) added. Circles and squares 
show objects with stellar- and gas-dynamical evidence for BHs, 
respectively. Upper limits on M, are plotted as crosses. The 
correlation may be only the upper envelope of a distribution 
that extends to smaller M,. 

An MDO detection in NGC 4486B is interesting because 
one promising explanation for double nuclei requires a BH. 
Double nuclei pose an interesting problem because they 
cannot reasonably be two star clusters in orbit around each 
other: the dynamical friction timescale for the decay of 
their relative orbits would be short. Also, for two close 
double nuclei to appear in the present sample of galaxies 
observed with HST, there should be many more at larger 
separations where merger timescales are still long. These 
are not seen. Tremaine (1995) therefore developed a model 
of M 31 in which both nuclei are parts of the same eccentric 
disk of stars. The brighter nucleus (PI), which is also 
farther from the baryecnter, results from the lingering of 
stars near apocenter, and the fainter nucleus (P2) results 
from the disk's increase in density toward the center. A BH 
near P2 is required to make the potential almost Keplerian. 
Only then can an eccentric disk be maintained by modest 
restoring forces such as the disk's self-gravity. 
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NGC 4486B now is a second example of a double nucleus. 
It, too, contains an MDO. This supports Tremaine's 
model. The details are different here; neither nucleus is 
at the galaxy center. If the double nuclei are an embedded 
disk, then the disk may have a different orientation or a 
central hole (Lauer et al. 1996). More generally, our results 
support models in which the survival of a double nucleus 
is connected with the presence of a BH. 
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